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Q ■ ABSTRACT 

We report the discovery of a dwarf Seyfert 1 active galactic nucleus (AGN) 
with a candidate intermediate-mass black hole hosted by the dwarf galaxy SDSS 
^ ■ J160531.84+174826.1 at z = 0.032. A broad component of the Ha line with 

FWHM=781km s _1 is detected in its optical spectrum, and a bright, point-like 
nucleus is evident from a HST imaging observation. Non-thermal X-ray emission 
is also detected from the nucleus. The black hole mass, as estimated from the 
luminosity and width of the broad Ha component, is about 7 x 1O 4 M . The host 
^ ■ galaxy appears to be a disk galaxy with a boxy bulge or nuclear bar; with an 

absolute magnitude of Mr = —17.8, it is among the least luminous host galaxies 
2 ' ever identified for a Seyfert 1. 



Subject headings: galaxies: active - - galaxies: dwarf — galaxies: individual 
(SDSS J160531. 84+174826.1) — galaxies: nuclei — galaxies: Seyfert 



1. Introduction 

Supermassive black holes (SMBHs) with masses Mbh ^ 1O 6 M have been convincingly 
inferred to be present in the centers of nearby galaxies (see Kormendy 2004 for a review), 
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and are generally believed to reside in most, if not all, galaxies with a spheroidal stellar 
component (Kormendy & Gebhardt 2001; Ferrarese & Ford 2005). However, little is known 
about their low mass counterparts, i.e. black holes with masses in the range of 10 3 — 1O 6 M , 
presumably in the centers of (dwarf) galaxies. These intermediate-mass black holes (IMBHs) 
may provide the 'missing link' in understanding the formation and evolution of SMBH seen 
today. In current models of galaxy evolution in a hierarchical cosmology, SMBH must have 
formed from much less massive "seed" black holes and grown up by accretion and/or merging. 
It is likely that there exists a population of IMBHs in the present universe, that have not 
the opportunity to be full-grown (e.g., Islam, Taylor & Silk 2004). 

The search for IMBH turns out to be a difficult task, however, since IMBHs are beyond 
the reach of direct measurement using star/gas dynamics, by which most nearby SMBH 
are unveiled (see van der Marel 2004 for a review). The most promising approach is to 
search for dwarf active galactic nuclei (AGN) that are hosted in small galaxies — a scaled- 
down version of typical Seyfert galaxies and quasars. So far, convincing evidence for the 
existence of IMBHs has been found in only two AGN, NGC 4395 (Filippenko k Ho 2003, 
hereafter FH03; M BH = (3.6 ± 1.1) x 1O 5 M , Peterson et al. 2005, hereafter P05) and Pox 
52 (M BH « 1.6 x 10 5 M o , Barth et al. 2004, hereafter B04), with NGC 4395 being the only 
one having accurate black hole mass measurement via reverberation mapping. In addition, 
Greene & Ho (2004; hereafter GH04) found 19 IMBH candidates from the Sloan Digital Sky 
Survey (SDSS; Stoughton et al. 2002) First Data Release (Abazajian et al. 2003), whose 
black hole masses were estimated by using the linewidth-luminosity-mass scaling relationship 
(Kaspi et al. 2000). All these objects, except NGC 4395, have a very high accretion rate 
close to the Eddington limit. 

In this paper, we report the discovery of a dwarf AGN with the estimated black hole mass 
as low as ~ 1O 5 M , that is hosted by SDSS J160531.84+174826.1 (hereafter J1605+1748)- 
a dwarf (most likely disk) galaxy at a redshift z = 0.032. It was found from our on-going 
program of a systematic search for AGN with candidates IMBHs or hosted in dwarf galaxies 
from the SDSS Fifth Data Release. We assume a cosmology with H =70 km s -1 Mpc -1 , 
f2M=0.3 and f2 A =0.7. 

2. Data Analysis 

2.1. Optical Spectrum 

J1605+1748 was spectroscopically observed in the SDSS on June 12, 2005 with 4048s 
exposure and was classified as a galaxy by the SDSS pipeline. It has a redshift of z — 0.03167 
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as determined using the [O III] A5007 line. Figure 1 shows the rest frame spectrum with the 
Galactic reddening corrected (E B _ V = 0.05, Schlegel et al. 1998), which is dominated by 
starlight of the host galaxy. To subtract the starlight and the nuclear continuum, we followed 
our method^ as described in detail in Zhou et al. (2006). The fit was good; the standard 
deviation of the distribution of the relative residuals f residual/ fsDSS in the emission line-free 
region around Ha and H/5 is ~ 0.05, just the noise level with respect to the signal. The 
left-over emission line spectrum is plotted in Figure 1. 

We fit the emission lines using the code described in detail in Dong et al. (2005). The 
peaks of the [N II]AA6548, 6583 doublet are well separated from the Ha line, that shows an 
apparent broad component, owing to the relatively high S/N ratio (> 20) in this spectral 
range and the narrowness of the lines. We model the [N II] doublet with two Gaussians 
with the line ratio A6583/A6548 fixed to the theoretical value 2.96. We use two Gaussians 
to model the narrow and broad components of Ha, assuming the narrow component has 
the same profile and redshift as the [N II] doublet (see Zhou et al. 2006 for a full account 
of this assumption). A good fit is achieved, with a minimum reduced \ 2 — 0-73 (d.o.f = 
92), yielding a line width 794 ± 42km s _1 FWHM for the broad Ha component (see Figure 
1). It has been noted by Veron-Cetty et al. (2001) that Lorentzian is a better profile to 
describe the broad lines in narrow line Seyfert Is, that have FWHMs < 2000km s -1 . In 
view of this argument, we use a Lorentzian for the broad Ha component and repeat the 
above fit; this yields a similarly good fit with reduced y 2 = 0.71 (d.o.f = 92), but a model 
over-predicted in the [N II]A6548 region by 7% of the [N II]A6548 peak height. So we adopt 
the Gaussian fit. To test the reliability of this broad Ha component, we model the whole 
Ha line with a single Gaussian with the center and width as free parameters; the result is 
unacceptable with a reduced \ 2 — 2.75 (d.o.f = 95) and large residuals remained around 
Ha. We therefore believe that the broad Ha line is real. For the H/3 line, since its S/N 
ratio is relatively low (~ 10), we fix the profile parameters to those of Ha, for both the 
narrow and broad components. For the [O III] AA4959, 5007 doublet, we only fit [O III] A5007 
alone with one single Gaussian because [O III] A4959 is jaggy. We find the fit is good and 
no extended wing for [O III] A5007. We fit other narrow lines with a single Gaussian profile. 
The results of emission line fitting are listed in Table 1. The Balmer decrement is found to 
be 4.2 for the broad lines and 4.0 for the narrow lines, indicating an extinction color excess 
Eb-v =0.36 and 0.31, respectively, assuming an intrinsic Balmer decrement of 3 (Dong et 
al. 2005, Zhou et al. 2006) and a SMC-like extinction curve. The ratios of the narrow lines 



1 Starlight is modeled with 6 spectral templates which were built up from the library of simple stellar 
populations of Bruzual & Chariot (2003) by using the Ensemble Learning Independent Component Analysis 
method (see Lu et al. 2006 for a detailed description). 
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[O III]A5007/H/3 > 3 and [N II]A6583/Ha > 0.6 place J1605+1748 into the AGN regime on 
the diagnostic diagram (Veilleux & Osterbrock 1987). It would be considered as a Seyfert 
1.8 in the Osterbrock (1981) classification scheme, just like NGC 4395 (Filippenko & Sargent 
1989, hereafter FS89) and POX 52 (B04). 



2.2. Optical Image 

We retrieved an archival HST image for J1605+1748, taken with the Wide Field Plan- 
etary Camera 2 (WFPC2) in June 1995 with a 500s exposure. The galaxy fell on the 
Wide Field camera (0".l/pixel) as a "bonus" in an observation of Mkn 298 (Malkan et al. 
1998). The F606W filter was used, which has a mean wavelength of 5947A and a FWHM 
of 1500A (1997 May WFPC2 SYNPHOT update). The data reduction was carried out 
following Malkan et al. (1998). On the HST image, galaxy is well resolved thanks to the 
superb spatial resolution. It is largely elongated with a major axis of about 6", and a bright, 
point-like nuclear source is evident. We perform structural decomposition using a "two-step" 
strategy: firstly fit to the one-dimensional surface brightness profile and then 2-D fitting to 
the image using the package GALFIT (Peng et al. 2002). We construct surface brightness 
models, for both the 1-D and 2-D cases, with a combination of different structural compo- 
nents: a point source represented by a point-spread function (PSF) that is generated with 
the Tiny Tim software (ver. 6.3)0 for WFPC2; galactic component(s) parameterized by ei- 
ther an exponential or a Sersic r 1//ra function (Sersic 1968, see Graham & Driver 2005 for a 
concise review), or a combination of any of these two. We start off the fitting procedure with 
the simplest one-component model, and add in one more component only if the fit can be 
improved significantly by doing so. An additional constraint comes from the above optical 
spectral decomposition, which sets an upper limit on the contribution from a point source 
to be ~ 15% the total flux at 6000A. 

Firstly, the azimuthally averaged 1-D radial profile of surface brightness is extracted 
using the IRAF ELLIPSE task. The best-fit model, which is plotted in Figure 2 along with 
the measured profile, is composed of a central point source plus two Sersic components, one 
dominating the inner part of the galaxy and the other dominating the outer part. Leaving 
off the point source or one of the two Sersic gives rise to unacceptable fits with significant 
and strongly structured residuals. The need for a central point source is not surprising, as 
revealed from the HST image (Figure 3). The presence of the second galactic component is 
found to be significant, as shown in Figure 2. Secondly, we perform 2-D image decomposition 
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using GALFIT, taking the results of the 1-D fits as initial values of the fitting parameters. 
The best-fit model turns out to be the same as for the 1-D fitting, but with somewhat 
different parameter values, especially for the central point source. We adopt the result from 
the 2-D fitting (Table 2) in this paper because 2-D modeling can recover a nuclear source 
much reliable than 1-D fitting (Peng et al. 2002). Also displayed in Figure 3 are the images 
of the residuals and of the model components. The fit is good in general, with the standard 
deviation of the distribution for the fractional residuals f residual / fusr ~ 0.1 within the sky 
region 10 times or more brighter than the sky level. In comparison, other models with the 
same or less number of model components give rise to unacceptable fitting results that have 
much larger \ 2 an d structures in the residuals. 

For easy comparisons with NGC 4395 and POX 52, We derive the Johnson magnitudes 
from the F606W flux using the IRAF/SYNPHOT package. We assume a continuum similar 
to that of NGC 4395 (/„ oc z/~ L5 , Filippenko, Ho, & Sargent 1993) for the central point 
source, a template bulge spectrum (Kinney et al. 1996) for the inner Sersic. We further 
make use of the SDSS spectrum (§2.1) for the outer Sersic since it accounts for > 90 per 
sent of the total F606W flux. The results are given in Table 2. 



2.3. X-ray Data 

As a remote member of the cluster Abell 2151, J1605+1748 was observed serendipitously 
several times by ROSAT and XMM-Newton. We retrieved from archives the observations 
with good effective exposures, one with the ROSAT Position Sensitive Proportional Counter 
(PSPC) in 1993, one with the ROSAT High Resolution Imager (HRI) in 1997, and one with 
the XMM-Newton EPIC-Pn| in 2003. X-ray data reduction were performed following the 
standard procedures, using the FTOOLS (ver. 6.1) utilities and the Science Analysis Sys- 
tem (SAS, ver 6.0.0) for ROSAT and XMM-Newton data, respectively. Table 3 lists some 
relevant parameters for the X-ray observations and data reduction. For the XMM-Newton 
data, exposure periods which suffered from high flaring background caused by soft protons 
are removed. Source counts are extracted from a circle (see Table 3) and background events 
are extracted from source-free regions using a concentric annulus for PSPC and HRI, and 
two circles at the same CCD read-out column as the source position for PN. In each of 
these observations an X-ray source was detected at the > 6a significance levels at a position 
coincident with the optical point source of the galaxy within one a positional error. The po- 
sition of the X-ray source detected with XMM-Newton is RA=16:05:31.87, Dec=17:48:26.11 



3 The source position fells into a gap between CCD chips on the two EPIC-MOS detectors. 
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(J2000), with an error circle of 0".9 (la) in radius. 

The spectral fitting is carried out using XSPEC (ver. 12.2.1). Due to the small number 
of source counts, 58 for PSPC (0.1-2.4 keV) and 79 for PN (0.3-8keV), the spectra are 
rebinned to have a minimum of 6 counts in each energy bin, and the Cash-statistic is adopted 
in the minimization instead of x 2 - We fit simultaneously the two spectra assuming the same 
spectral models with the parameters tied together except the normalizations. The spectra 
can be fitted with a power law with a photon index T = 2.1 ± 0.3 and an absorption column 
density close to the Galactic value (3.5 x 10 20 cm~ 2 ). The unabsorbed 2-10 keV flux is 
3.0 x 10~ 14 ergs s^ 1 cm -2 measured by XMM-Newton, corresponding to a luminosity of 
7.0 x 10 40 ergs s _1 . The best fit normalization at 1 keV is a factor of 2 higher for ROSAT 
than for XMM-Newton. To test the significance of this variation, we search for the maximum 
confidence contours in the parameter space of the two normalizations (two free parameters) 
within which they differ from each other; we find that the variation is marginally significant 
at the ~90% level. The HRI flux is found to be roughly consistent with the XMM-Newton 
value assuming the same spectral shape. 



3. Results and Discussion 

In summary, observations of J1605+1748 reveal the presence of an optical bright point 
source based on its HST image, a broad Ha component with FWHM 781 km s _1 (after 
correction for the instrumental broadening of 141 km s^ 1 FWHM) and AGN-like narrow 
emission line ratios, a non-thermal X-ray luminosity (7.0 x 10 40 ergs s _1 in 2-10 keV) and 
possible X-ray variability. All these observational facts point to the presence of a Seyfert 
1 type dwarf active nucleus in J1605+1748. This can be further supported by the striking 
similarities between J1605+1748 and the two well known dwarf AGN NGC 4395 and POX 
52, whose properties are summarized in Table 4. In fact, the AGN in J1605+1748 lies in 
between NGC 4395 and POX 52 in terms of the observed luminosities of both Ha and optical 
continuum. Furthermore, J1605+1748 is more luminous at hard X-rays than NGC 4395 (no 
data available for POX 52 so far). The equivalent width of its broad R(3 component (~ 20A, 
relative to the nuclear continuum calculated from the HST F606W flux of the point source) 
is typical of Seyfert 1.8/1.9 nuclei and similar to that of NGC 4395 (FS89). 

We estimate the mass of its central black hole following the widely used linewidth- 
luminosity-mass scaling relation (e.g., Vestergaard 2002; McLure & Jarvis 2002; GH04). 
The intrinsic monochromatic luminosity at 5100A is estimated to be L5100 = ALa(5100A) = 
4.0 x 10 41 ergs s _1 from the Ha luminosity using the L Ha —L 5100 relation given by Greene & Ho 
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(2005B incorporating the extinction correction that is derived from the Balmer decrement 
for the broad lines (§2.1). Using Equation (6) in Greene & Ho (2005), we find a black 
hole mass 6.9 x 10 4 M Q . The uncertainty of the black hole mass thus estimated is not well 
understood, however. As pointed out by Vestergaard & Peterson (2006), the statistical 
accuracy of the masses from such a scaling relationship is a factor of ~ 4 (la), and, for 
individual mass estimates, the uncertainty can be as large as an order-of-magnitude. Using 
the BLR-size— luminosity relations from Kaspi et al. (2005) and Bentz et al. (2006), we find 
a black hole mass of 5.8 x 10 4 , 2.6 x 10 5 M Q , respectively; here we adopt the formalism and 
scale factor from Peterson et al. (2004) and Onken et al. (2004), with L5100 = 4.0 x 10 41 ergs 
s \ ciinc — 332 km s _1 estimated from the broad Ha component. 

Assuming a spectral energy distribution (SED) similar to that of NGC 4395, the bolo- 
metric luminosity for J1605+1748 is estimated to be ~ 3.7 x 10 42 ergs s" 1 by scaling 
their L5100 luminosities using the NGC 4395 data in P05. This is consistent with the bolo- 
metric correction L^oi ~ 9L 5 i o used for normal QSOs (Kaspi et al. 2000). Then, the 
Eddington ratio L^/L^d is ~ 0.3 for a black hole mass of 6.9 x 10 4 M Q . This indicates that 
J1605+1748 is at a high accretion state, similar to POX 52. The broad-band spectral index 
a ox = —0.3838 log[F !/ (2keV)/F l/ (2500 A)] ~ —1.1 agrees well with the extrapolation down to 
low luminosity of the a ox — luminosity relation given by Strateva et al. (2005, see also Yuan 
Brinkmann & Siebert 1998) for radio-quiet AGNs, similar to the result for a small sample 
of IMBH AGN (Greene & Ho 2006). 

As results from the above HST image decomposition analysis (§2.2), the galaxy is com- 
posed of three components, a central point source as the AGN, an inner Sersic with r e = 0".18 
(114pc), and an outer Sersic with r e = 1" .58 (lkpc), which account for 6.0, 2.6, 91.4 per 
cent of the total F606W light, respectively. We tend to identify the outer Sersic to be a 
galactic disk, because of 1) its exponential-like radial profile with the Sersic index n = 0.8, 
close to 1 for exponential, 2) the disky isophote shape (c = —0.1), 3) the apparent extreme 
flattening (Hubble type E6.3 based on the axis ratio), that is very unlikely for an elliptical 
(Binney & de Vaucouleurs 1981), and particularly, 4) the existence of the additional inner 
Sersic component. 

With an effective radius of 0".18, the inner Sersic is merely resolved by the Wide Field 
Camera; however, fits without this component result in large residuals and an over-predicted 



4 This value is well consistent with that estimated using the HST F606W flux of the point source assuming 
a power-law spectrum (fx oc A~" A ), (3.5 — 4.6) x 10 41 ergs s _1 for a\ = — 2. 

5 This somewhat shallower inner profile may be caused by dust extinction in the disk, which is indicated by 
the Balmer decrements, since the extinction aggravates inwards for an inclined disk (i = cos -1 (b / a) = 66°). 
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point source, as displayed in Figure 2. Its size is too large for being a nuclear star cluster. 
The fitted image has a boxy shape (c = 0.6) and a rather flat profile (n = 0.4), suggesting 
a box-/peanut-shaped bulge or an edge-on nuclear bar — their presence in disk galaxies are 
known for many years (eg. Burbidge & Burbidge 1959; Evans 1951; de Vaucouleurs 1975). 
We prefer the nuclear bar speculation because there is a ~17° mis-alignment between its 
major axis and that of the disk, that is more natural for the bar scenario; moreover, a bar 
can provide an efficient mechanism of fueling gas to the AGN but a bulge cannot (e.g., 
Shlosman et al. 1989). In fact, a boxy bulge may evolve from a bar, or may simply be (part 
of) a bar viewed edge-on (see Kormendy & Kennicutt 2004 for a review). Regardless of its 
nature, we predict the mass of the central black hole based on the empirical Mbh ~ ^bui 
correlations using the absolute magnitudes (Table 3, corrected for the internal extinction 
Eb-v — 0.36). Using the formalism from Marconi & Hunt (2003) and McLure & Dunlop 
(2002), we find 2.5 x 10 5 and 2.6 x 1O 4 M , respectively, roughly consistent with our estimate 
from the scaling relation considering the large uncertainties. Furthermore, following the 
empirical relation between Mbh and bulge concentration (Graham et al. 2001, Graham et 
al. 2003), we obtain Mbh — 1-6 x 1O 5 M , also consistent with the virial mass estimate. 
Either bars or boxy bulg kind of "pseudobulges" [Kormendy 1993]) form out of disk 

material (Kormendy & Kennicutt 2004). So, probably just as stated by Carollo (2004), "the 
requirement for making a supermassive central black hole is that the galaxy is capable of 
reaching sufficiently high central baryonic densities." 
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Table 1. Fitted emission line parameters 



Line 


Centroid a 




FWHM b 


Flux c 




[0 II]A3727 


3728.47 ± 0. 


.51 


281 ± 102 


28.0 


± 


8.5 


H/9 (narrow) d 


4862.14 ± 0. 


.24 


137 


17.0 


± 


2.9 


H/3(broad) c 


4860.93 ± 1 


.26 


781 


27.3 


± 


5.8 


[0 J//]A5007 


5008.25 ± 


.08 


160 ± 12 


53.7 


± 


2.8 


[0 7]A6301 


6302.83 ± 


.86 


277 ± 84 


9.3 


± 


2.6 


Ha (narrow) 


6564.53 ± 0. 


,06 


137 ± 8 


67.3 


± 


4.1 


Ha (broad) 


6565.78 ± 


.30 


781 ± 42 


115.2 


± 


5.3 


[iV J/]A6583 d 


6585.22 




137 


38.4 


± 


1.9 


[5 7/JA6716 


6718.43 ± 


.18 


142 ± 20 


18.2 


± 


1.9 


[5 //]A6731 f 


6732.81 




142 


14.6 


± 


1.9 



Note. - - (a): Vacuum rest-frame wavelengths, in unit of 
A. (b): Corrected for the instrumental broadening; in unit 
ofkms -1 . (c): In unit of 10~ 17 ergs s~ x cm~ 2 . (d): Adopting 
the profile of narrow Ha. (e): Adopting the profile of broad 
Ha. (f): Adopting the profile of [S II]A6716. 
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Table 2. Two-Dimensional Image Fitting Parameters 



Func. 


m 606 a 






r e 


n 


b/a 


PA 


c 


notes 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


PSF 


21.1 


-14.6 


-14.6 












AGN 


Sersic 


21.9 


-13.9 


-12.4 


0.18 


0.4 


0.28 


82 


+0.6 


Bar? 


Sersic 


18.1 


-17.8 b 


-16.4 


1.58 


0.8 


0.37 


65 


-0.1 


Disk 



Note. - Col. (1): Components used in the fit. Col. (2): The 
F606W integrated magnitudes on the Vega system. Col. (3): The ab- 
solute Johnson R magnitudes. Col. (4): The absolute Johnson B mag- 
nitudes. Col. (5): The effective radius of the Sersic law, in arcsecond. 
Col. (6): The Sersic exponent. Col. (7): Axis ratio. Col. (8): Position 
angle, in degree. Col. (9): Diskiness (negative) /boxiness (positive) pa- 
rameter, defined in equation (3) of Peng et al. (2002). (a): All the 
magnitudes are corrected for the Galactic extinction; see the text for 
the transformation from mag606 to Johnson B and V magnitudes, (b): 
In good agreement with the SDSS r— band absolute magnitude derived 
from the fit with an exponential model, M r = —17.8. 
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Table 3. X-ray observations and data reduction information 





XMM/EP1C-PN 


ROSAT/HRI 


ROSAT/PSPC 


Obs. id 


0147210301 


RH703861A01 


RP800517N00 


good exposure 


5.7ks 


54.9ks 


12. Iks 


-^source extraction 


30" 


25" 


45" 


net count rate* 


1.2 ±0.2 


0.12 ±0.02 


0.35 ±0.08 


flux* 


3.4 ±0.6 


4.2 ±0.7 


7.1 ±1.4 



Note. - (*): Measured in the 0.3-2.4keV range; count rate in 
unit of 10~ 2 cts s _1 , and flux (without Galactic absorption correction), 
10~ 14 ergs s' 1 cm' 2 . 
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Table 4. Comparison with NGC 4395 and POX 52 



Object 


M Host 


M^ GN 


j br 


Lx,Soft 


Lx,Hard 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


NGC 4395 


-17.5 a 


-10.8 a 


1.2 b 


1.3 C 


1.5 d 


J1605+1748 


-16.4 


-14.6 


27 


1100 c 


7.0 


POX 52 


-16.8 f 


-16.7 f 


120 f 







Note. — Col. (4, 5): In unit of 10 38 ergs s _1 . Col. (6): 
Unabsorbed X-ray luminosity in the 2 — lOkeV range, in 
unit of 10 40 ergs s -1 . (a): From FH03. (b): From FS89. 
(c): Observed X-ray luminosity in the 0.2 — 2keV range 
from Moran et al. 1999. (d): From Vaughan et al. 2005. 

(e) : Mean observed X-ray luminosity in the 0.3 — 2.4keV. 

(f) : Form B04. 
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Fig. 1. — SDSS spectrum and the line fit. Left: The SDSS spectrum (top) and the 
starlight/continuum subtracted residual (bottom, offset downward by 2 units for clarity). 
Right: The line decomposition in the Ha+[N II] region. The upper panel shows the original 
data (thin solid line), the fitted narrow lines (dotted lines), the fitted broad Ha (thick solid 
line) and the sum of the fit (gray dashed line). The lower panel, the residual of the fit. 
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Fig. 2. — One dimensional fitting of the radial profile, with PSF + two Sersic model (left 
panel) and PSF + a Sersic model (right panel). Each panel shows the point-source component 
(dotted), the host-galaxy component (s) (dashed), and the sum (solid). Residuals are also 
displayed in the bottom. As indicated in the right panel, fit with only one Sersic for the 
host galaxy yields large residuals left and an over-predicted point source component. The 
results of the 1-D fits are taken as initial values of the fitting parameters of the subsequent 
2-D decomposition. 
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Fig. 3.— Two-dimensional decomposition of the HST/WFPC2 F606W image. Top: the 
original data (left) and the residual from the best-fitting PSF + two Sersic model (right). 
The standard deviation of the distribution for the fractional residuals f residual I I'hst ~ 0.1 
within the sky region 10 times or more brighter than the sky level. Bottom: the best-fit 
inner Sersic component (left) and outer Sersic component with the contours of the original 
image overlaid (right). Images have the same size. 
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Supplementary information: 

[See more at http://staff.ustc.edu.cn/~xbdong/ J1605+1748/] 



Dwarf Disk Gaiaxy with an iMBH 



J160531.8AH74826.1 : B=-f6.4. R=-f.7,3 
z=0.032 , ■ 



4 A 
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cluster (Abel! 21S1) 
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Fig. 4.— HST/WFPC2 F606W preview image for J1605+1748 and its neighbor Mkn298. 
They are in the Hercules cluster (Abell2151). Displayed with a logarithmic stretch. 



